(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     33608,        637]
NotebookOptionsPosition[     33034,        615]
NotebookOutlinePosition[     33381,        630]
CellTagsIndexPosition[     33338,        627]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"data", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "1.91373"}], ",", 
       RowBox[{"-", "0.799904"}], ",", "2.04001"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "0.935453"}], ",", 
       RowBox[{"-", "0.493735"}], ",", "0.959304"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.630964", ",", 
       RowBox[{"-", "0.653075"}], ",", "0.728477"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.310857", ",", 
       RowBox[{"-", "0.018258"}], ",", "0.301885"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0431084", ",", "1.25321", ",", "1.19012"}], "}"}]}], "}"}]}],
   ";"}]], "Input",
 CellChangeTimes->{{3.609308870343686*^9, 3.609308961165881*^9}, {
  3.60960770376698*^9, 3.6096077329416485`*^9}, {3.6096078143023024`*^9, 
  3.6096078847353306`*^9}, {3.609778220829998*^9, 3.6097782230671263`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "  ", 
   RowBox[{
    RowBox[{"Non", "-", 
     RowBox[{"Weighted", " ", "Fit", "  ", "fit"}]}], "=", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"data", ",", 
      RowBox[{"Norm", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", "y"}], "}"}], "-", 
        RowBox[{"{", 
         RowBox[{"x0", ",", "y0"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x0", ",", "y0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "y"}], "}"}]}], "]"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"fit", "=", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"data", ",", 
      RowBox[{"Norm", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", "y"}], "}"}], "-", 
        RowBox[{"{", 
         RowBox[{"x0", ",", "y0"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x0", ",", "y0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "y"}], "}"}], ",", " ", 
      RowBox[{"Weights", "\[Rule]", 
       RowBox[{"1", "/", 
        RowBox[{
         RowBox[{"data", "[", 
          RowBox[{"[", 
           RowBox[{"All", ",", "3"}], "]"}], "]"}], "^", "2"}]}]}]}], "]"}]}],
    ";"}]}]], "Input",
 CellChangeTimes->{{3.6093090050093884`*^9, 3.6093090050153885`*^9}, {
  3.609309130355558*^9, 3.6093091338597584`*^9}, {3.609309248294304*^9, 
  3.6093092895736647`*^9}, {3.609589968883603*^9, 3.609589992731967*^9}, {
  3.6096075061266756`*^9, 3.609607506782713*^9}, {3.6096076009530993`*^9, 
  3.6096076418674393`*^9}, {3.6096077565059967`*^9, 3.609607810889107*^9}, {
  3.609778108552576*^9, 3.6097781390073185`*^9}, {3.609778178536579*^9, 
  3.6097782063671713`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ellipsoid", "=", 
    RowBox[{"fit", "[", 
     RowBox[{"\"\<ParameterConfidenceRegion\>\"", ",", 
      RowBox[{"ConfidenceLevel", "\[Rule]", "0.95"}]}], "]"}]}], ";"}], "\n", 
  RowBox[{"(*", 
   RowBox[{"fit", "[", 
    RowBox[{"\"\<ParameterConfidenceIntervalTable\>\"", ",", 
     RowBox[{"ConfidenceLevel", "\[Rule]", "0.95"}]}], "]"}], "*)"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.6095900263808913`*^9, 3.6095900822370863`*^9}, {
  3.609777656365713*^9, 3.6097776612219906`*^9}, {3.609777747917949*^9, 
  3.609777794820632*^9}, {3.6097778625935087`*^9, 3.609777883508704*^9}, {
  3.609778074478627*^9, 3.609778077062775*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"source", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"solution", "  ", "=", " ", 
   RowBox[{
    RowBox[{
    "fit", "[", "\"\<ParameterConfidenceIntervalTableEntries\>\"", "]"}], "[", 
    RowBox[{"[", 
     RowBox[{"All", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Graphics", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Opacity", "[", "0.2", "]"}], ",", 
     RowBox[{"EdgeForm", "[", 
      RowBox[{"Opacity", "[", "0.75", "]"}], "]"}], ",", "White", ",", 
     RowBox[{
      RowBox[{
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"Most", "[", "#", "]"}], ",", 
         RowBox[{"Last", "[", "#", "]"}]}], "]"}], "&"}], "/@", "data"}], ",", 
     RowBox[{"Opacity", "[", "1", "]"}], ",", "Red", ",", "ellipsoid", ",", 
     RowBox[{"PointSize", "[", "0.0125", "]"}], ",", "Blue", ",", 
     RowBox[{"Point", "[", "source", "]"}], ",", "Red", ",", 
     RowBox[{"Point", "[", "solution", "]"}], ",", 
     RowBox[{"PointSize", "[", "0.0083", "]"}], ",", "White", ",", 
     RowBox[{
      RowBox[{
       RowBox[{"Point", "@", 
        RowBox[{"Most", "[", "#", "]"}]}], " ", "&"}], "/@", "data"}]}], 
    "}"}], ",", 
   RowBox[{"Background", "\[Rule]", "Black"}], ",", 
   RowBox[{"ImageSize", "\[Rule]", "600"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.609590106915498*^9, 3.6095901069304986`*^9}, {
   3.6095901404994183`*^9, 3.6095901822288055`*^9}, {3.6096045462793818`*^9, 
   3.609604554893875*^9}, {3.6096051938874235`*^9, 3.6096051981026645`*^9}, {
   3.6096052736959877`*^9, 3.609605289190874*^9}, {3.609605930217539*^9, 
   3.6096059408731484`*^9}, {3.6096074563878307`*^9, 
   3.6096074593860025`*^9}, {3.609607960146644*^9, 3.609608000193935*^9}, {
   3.609777637822652*^9, 3.609777676526866*^9}, {3.609777888303979*^9, 
   3.609777895160371*^9}, 3.609777953114686*^9, {3.6097780487801576`*^9, 
   3.609778050608262*^9}}],

Cell[BoxData[
 GraphicsBox[
  {GrayLevel[1], Opacity[0.2], EdgeForm[Opacity[
   0.75]], {DiskBox[{-1.91373, -0.799904}, 2.04001], 
    DiskBox[{-0.935453, -0.493735}, 0.959304], 
    DiskBox[{0.630964, -0.653075}, 0.728477], 
    DiskBox[{0.310857, -0.018258}, 0.301885], 
    DiskBox[{0.0431084, 1.25321}, 1.19012]}, 
   {RGBColor[1, 0, 0], Opacity[1], 
    InterpretationBox[
     GeometricTransformationBox[
      CircleBox[{0.03003887748055729, -0.12245764066673417`}, \
{0.3552833933340725, 0.14207314729546827`}], {{{-0.30031693031193113`, \
-0.9538394735845328}, {0.9538394735845328, -0.30031693031193113`}}, Center}],
     FittedModels`ParameterEllipsoid[{
      0.03003887748055729, -0.12245764066673417`}, {0.3552833933340725, 
      0.14207314729546827`}, {{-0.30031693031193113`, 
      0.9538394735845328}, {-0.9538394735845328, -0.30031693031193113`}}]], 
    {RGBColor[0, 0, 1], PointSize[0.0125], PointBox[{0, 0}], 
     {RGBColor[1, 0, 0], 
      PointBox[{0.03003887748055729, -0.12245764066673417`}]}, 
     {GrayLevel[1], PointSize[0.0083], PointBox[{-1.91373, -0.799904}], 
      PointBox[{-0.935453, -0.493735}], PointBox[{0.630964, -0.653075}], 
      PointBox[{0.310857, -0.018258}], PointBox[{0.0431084, 1.25321}]}}}},
  Background->GrayLevel[0],
  ImageSize->600]], "Output",
 CellChangeTimes->{
  3.609778496364758*^9},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJzVnXm4lcWR8N/cc+69rKK4ISqi4oYbbkSNC2JcQFGMKy4ZCIogwQ3cl4BG
jBr3gEFQBMXBDRRRFHeNUaMmMclnYpxEzWicSTIZJ8lkknyZzNe/quq3u+7p
i3m+/+Y8nuPlnLe6lq6urq6q7h4zYcbk06ZNmDFl4oTBh5w74ezJUyZOH3zw
WeeGrxqfqaq2E8L7mMEVf/9PVcWP9OeJVZW+/t/wxY/lszkkfAyz9z7h3RjP
140b2tvb53R0dPCu9LN9dWdn57ez9/PhXfXQ33geyPgeTittw6StPlX2EhJO
0r/flM/OPoZ+pMFe3Gg0aHU+LTfAVDXvW2utte5aZ511eN+23nrrNWYOHhx6
44vbbrvtuO2224531SGf7fvusssuu2WvPXfdddeqh/7G80DG94J11123aqPR
qvl0z549I2vLwrtNOZv0mc98JnIElR01F0crF9+Wz7UG2DM8CwzQi8Kb5mga
JLdusMEGoJ221VZbQdCx22+//VBorPodvNNOO0VOLtl8880jdYA9b1TNtw7J
SRpUk3O4kvO8fPbpMHKPCO9pbW3CDO2s6N27N+3SPnjAN3aHHXb4nJDRB4q6
kvBQnz59QE8TM0LfHGNN94mYD1bMT8ln7z6mSVCIWiwzVYGT6wYOHHj+Flts
AQrpl0MFbS/+fc6WW2559SabbMJzCGyloRxvqtE/ojtA0a2Sz55DjEm0ZqnJ
Gg2BerobwVZHC5qeB+y8886g4jf6YlWvXr0gb1az2YSrobUw91UcKxRHFGZk
iZGARsIOfXnYjjvuWJ2iOMA3eeutt4aVe/r16/d8xsZwk5ri2FNxLJfPHnTk
MdZ+FFfkAZU+TZrvgdQiC3QO7KIXjOohseU9tOUHdHx1WIfQd3BLf0IdLVdn
SbOd9PuVgwYNgisEjzRHOmp30TaXaptdqUXidKBI+zxtE2kjHSTdlUxtcydt
c4l8dgwxVe1K42XSXEdOIo/waM3x9trS4taWeByVg7IrtSX+5KtiS9tqSwvX
3NI15ZYmhJE5KLa0lba0QFvCPpxeaOnG1NJXNttsMxhfai0NiC1toS3dLp/t
I03uDOWLgxqIzOdIM+1oIs3wG8+MrGU9WNu4TduI1NDX9I/ownxtg7/5jt94
ZkDdxqbaxpwyHfR3tVDZ4W++y+noiM0M1GZu0WYY11/qSsoST8pjPXr0ODOo
ZCJlgLZxk7aBzFHZR8KAvjzoiJByn7bB4DkvmDMhMzyT1G8DbeMGbWOYsUMP
oLZin5ZJG03+TacxkIdH8PUU/Ovy2Tw0fGDvGbeQ/ahCMlYZp/TpoRGyv0Je
o4gZYyeF7kYANw4YMABzXK1K/cl3/MYzaTz20zauTm0gxIfDOLti0003FQE8
pW3wN9/xG8+kNvpqG1dpG/TPydYGQkRw1XNJiHzHbzyTZsHe2saVqY1jA53L
Q4fNCuNX2ngptcF3/MYzqY2e2sas1jZmb7zxxsLLK4kXvmtto1PbmKm90RX+
de2NruAC2aGQlyskakBnRcjveUh+qjWgqZCXKiQq+bWgQGgA/fZDheRPNADd
qm1Dm0JerJBQMDmo5oNBRSdus802P1ZI1AhtRPPrybaqLlIghhR8XBSGO7r6
UwFqoOr00pGBvepvCnGBQqD5GGFGFPz8TNHQL9dutNFGaC/c/1WBzlMgEEMA
6o8J/kCBUFIYvSRMnajzXxRohnw2kND9gQiIqT5SwqaG6eGBMD5HQdif9fFz
FUd/kxxigrCPE2E3bLjhhkiOR/5Lgc5RHLRDewis+rXiQGLYAFS0+qM+fpbi
GGTM0yAN/5viQHIoJT/xyB8U6EzFQTtL+/btC+PVfygEf/MdHUbH/U4hvuzR
0CZt/z6hoaMWBPXhkf9QoKmJlUWhowXNfykrCBg0x8HKJ/r4lNQlyD2y8ueE
46rgkURx/VaBJisOfJt7QpfQDdV/J3EhQvHW/00fP10fnxi4w32RyeV/Un/Q
ReCG818rxGkKMcA4vyyYGoVopo5fHPgW7+pXCnNq0pN7g9IL4wLToTAMmjsD
+rFw/68KMzFRBhViWgSmR+L/mmArrwj8Q97HCjZBwfqYjiXy+iiqSUOGDKFb
ZEb5pcKMV5ghJrMJwfwrzFqpc6BbFjgfKcw/KMyBoS/vDmMlCW5thYFeRozo
5ocKc4p2KNSeF4R6y/rrrw8XAraugqFI0C1W+xcKdnJiaW6g/CxGmcBsoDDQ
e2egW1j6QGF0zdUYZCwliQ9IEocl6aX3FUaXig24RD0TSwMTS/QSPFfvKcw4
hRkduLw9cJt6aROFgT/4PCV0Y/VzhTlBYfjq+mCbZOYVmM0UZmYwvMhHHvsn
hTlOYRDdpUE8PKIwWygMQrkxqKPI7V2FOTbJjZ/OCC6YwgxRmFOCKbktmFSR
208V5hiF4St+knWkwGydVOH2IDdxK95RmKMTDD/JhC4w2yoMgrwj6IjA/ERh
vpD6B5gk66FJ1jWMLt6ro5Ka+v7ZIcGgi3uicm8rzNgEQ3Opf3ZK/TOnf//+
AvN/FOZIheErDzPMw9Dp1Y8U5ogEw08JZtcEQ2cLzA8VZozC7GcwSQ/2UJj6
eV5vKczhCQ/6lmA+W4D5vsIc5mlLMHsVYL6nMKMVhrHg8ejKuYFOyzjg9V2F
GdWdrPcpyNpgDu0OZr8CjEZQqkOSvjGEE8yIBHN31J03FOZgr6NJd0YW9O11
hTmoO5jPF2C+ozCfVxjGHOMnjYWD01igMRlzrzkYxihDUh4UmEPSmMNki+f2
qsIcqDCYCYZ+ghmVxjYNymOvKMxIhcHuYD8SzGHJhtR2R+M8Fn0QM0+XJ5gx
yVaJbeP1ssKMUBhMK+Yy6c6RCjM3TGf7RX37lsLs312fjk19SmNDM5j9kqyx
8al/vpD6594o65ccDCsbbFvi52iFYV5YHGX9osJodKQxJYjmwrAoSDDHKgzz
D3OWPPaCwuyjMLAJuwnmeIXBczkuyuB5hfmcwmATsW9JBuMUhul0WGWv5/R/
eysM5EK2hGYE5sSCDDwMJEN6ou0khWGynBtl8KzC7KUwkAzpCUbjQOL1XBJ1
5xmF0YhPwYb8g8LgvYyNMnhaYT6rMDh0S4KrJoshgRlf0AMP08rPhEKfathO
I7ZVg2UDi4EEM7HAz2qF2SP1D65q4ue0gl4/qTC7p/5hCCd7MCn1Dz5p/1YY
liEzBg8enGg7XWFQ3DlxbD+hMLslWaMmCWZKQd80lljt2p2sz0iyru2oh2mV
9dSCrB9XGI2liduPD59gpikMK4XJ0YY8pjDDuqNNI3iyPhtS2Wulg8EhxstN
eAwG4S+KcjOYnRUGdVjgxuk5CoM/enKU26MKs1PqU1akibZzk9wejGPOwxCU
PDv4YQnPdIVhbXJDtCEah612VBhZygSVTzAacWwgzFFxDn5EYXZItHl7cL7C
QOySqG8ehpUGy4WE54LCWHhYYTQCKfYaVtNYuEhh7oljlNdyB4PuTHTz3EWF
/jEYbaOg15coDIP3yCiDZQ6GpmiyBQYZsFzvk8FspzBESWa7uVHDsRJnODPq
6EMKo7FTmX9oLunB5YXx42Fo6jw3b3+loAcPKsw2CtNqd2YpDNHofaKsH3Aw
rbI2GCbL+4LydmQwWysMTdFkgtFAcgPBHBtlfb+DaZX1lUnWyyMeg9mqOzxX
FfDc52Ba8VxVwGMwQ7rDc3UBz1IH04rn6tSnj8TxYzBbdofnmoLu/KODacVz
TQGPwWzRHZ7rCnjudTCteK4r4DGYzbvDc73CEOk+PdrEJQ6GiAjuQwuM6x+D
GawwiBxznmBuLPTPPQ6mVa9vLNgQg9lMYbBVRMUTzM2F8XO3g2kdpwbj5mCD
GaQwiBKRpnF6q8K4OWuxgyHYfYWzvQZDbm1WtL0Go3mJwvzzjUKfehiWqdc5
/21OoU8XOZhW3ZlT6FOD2aS7Pr2tIOu7HEyrrA3G2USD2TjJgGxJ8pXnKcyK
OP/yWuhgWHKSZEh45hVk4GFa1wvzCjIwGE34iL0mnp1g5hdkcKeDaZXB/IIM
DGaj7vTgjoIM7nAwLAWJiyc8dxRk4GFQQ9SxBYa5hDxfDjMg0fZUoC2NBc2Y
ef9ggYNppW1hgTYP00rbwgJtBrNhd7QtKtA238G00rZjGnMoCKmzyso15oTF
4v3BhvGu2uSz/dWwiH6rf//+8f0GJRlN+a3z9mAA8nINMsy0Kco9riADT1ur
DNQ89mTwxlKQm4KvDTboeDPo4aPBUV7Cwrkxb6+99gr+1MiRIy857LDDeFdt
8tl+6gknnHBy9hp/4oknVk35re/sz3/+80DG9+JAK22+FoQEb+BCDuDGYEsg
9uhC3xgvG3TXN1sKjAQGYBI+ngyDGBmuDn8vD9+B/4YRI0ZAGTRXG0Nt1Rt6
+Q7eeObBINNVwQkH9tngQ0Mf/SYD6MiCDtzuaGvVAY3bCrHweG2gZ2XwS5HB
0sAvNJ111FFHbSfk9Jh69NFHQwpifzmwABuwA1uwd1ihpz0FrT2t0eZOtOWq
ALc8mHh6F05v2W+//aYfeeSR1S6CvhNKoAjKoBBKoVgSS4cWesZwr596xltc
nUwkpT8rmC4kigbM3WeffS4YM2bMXoK2AxIgBZJeCmjhmNWJxcUwlIOivOc5
jK2+vy4hOyCBMQK3LwY9++aee+6JcPdXjJOPO+44GH0mLOGQMSUIdPGBipF6
i2PKGBH7xW7lrYEvcTkYpd8KPKLp8FcdpOiQ68Ldd9+dLl0S2o5JdQtKOodB
V9GStoI7yBhfaclK1VxD+Zn9Rl6f5080eDGCYx0LrZN2/yQ0uok+oC8YwAyQ
0cJF+1dGjRqFYmIfrotVDlks9qk41+miUcqfYiUbajnHasCoNaLTkHK0ZVoe
0E3Jmv7WG6nxPKUvwDOHxWKwr1nVXSyPsiC1VSYU/DtVGcmiw/LzQdNv++xn
P4tiVEcqv7B+U9BUeg07IEHEfZVfN0/rgqpnLPeiSgYZU+kT69wYksKFhpMk
6YTG5PVXJNRnWLWR8HCw46HVP+uXeKA/Xgg/woPYt2OEhybG5EWzswwpC9Ej
8hlxOaB+Zyf/o55uvpEOi4wLiUupiZWYAbCRZNiEXQmvHOjIbXXD+yq5iBFy
ECvilSnjhEQuY/KuMJDAu1fB5uhyQqqIcI2XWhURdXpIVdcXouBQj6Iss/Id
jMgIJXJdJbJ1LlGXrBmJRHYQdYrSd8XBBx/MlBLp27NgodQ1lmwbXYY4UT1c
eMS3SzKJ2JFY5XSMdfq+jr5W+9ZL6ZNkfeAYYUHfBE8f6sBQGl6wZ+q2St0T
xoKuRq1k9tVMkdTaAEMX048y233OUUavn++iVT2VMryJx0wT6dpThbIGPf1U
+JqBUO1RsHvqSUplxteslIoBr+OtiSYx786xyoy9HDGtFq1TwRga9BUuAEb4
dBUTo5w5Dg8EXncrGDH1nJpDjR5UDCIO0IYhnY6hdxnzGme2ciYJ+yDXrvQ0
Lg0SY5jST9VUJQb3g26cHboTCiwr6CL66yVikACSgJiDlBi6jtG6yGqy9nDE
tMYv2xMYk9/te+yxBz11pvbUN/bee+/nwtdSfbZLQYG03bocBmIQh2aXpBwm
Fvnw5G6OGIr5qMpLxFhpxBVBMjhe9FJ1Tuomeu7moM1IZmd9FM7rnMo62i42
sy6zsfyTK4HZ1ZHRumZrUxiJ/wWtQHur81QgXw+uCRotAVPLDTO4VkeLpAVo
UvAC17Go6MhsqAb5LbXCj2GOEhqJk55S8pk0OFYFMPGlL1BKEM6TwSqdySDa
qWB7tIytLtaRLjxayYB6uopoLV1nFacmv9ZF+WdSv+CyycR4kfYLynJPaA4D
CXeWY2fhfGikRIvhpJ5pTlYEdJw+StiYsgyZenZ0lKBmqJs8+LdCr1ymsrjm
gAMOoFdkelKj38BgXByntV5J/cmiyOxrKyVkQfRVMkw7OPSAT4vZi/8udMVX
Ulc8FvRyOos6K0xwASat5GtMzmt0TkoaEctz6IqhjgJmK2YHefav2jIdvixg
kl6YpRTQC3cFwVyJ+2QUPB1D17yUDok5YaZIpVSWXiMHQl5UMGt8u1pbH3eu
kqGHSTpcpuqvKvq5w4cPnxtHslVzOK+oPWkW8pf5xbJusEemTjJi2zr0LrL4
f/VxOIRTQTxb0X8teO/zTf+qbey58O8TI3oNi/lKjYnJzpOQk8KGbRx6wK+M
idy/6ONwCbeC+GuK/qoDDzxwUawMtfoXF49THZTkeF30YXk9airqvPhWDr1z
mv6UqIVTOK6uU/SXHXLIIQ+GeU2a3Sop9dJoknRp2EDBYylVNVmfIy0/JZI5
xKFfamkhefC/kqWDUziurlf0LNgeD82IR6zopfbd6kerU0PzKHhlmT9S45dG
i63hbTOa4oEB2hUl3MFldaOiZPG2PJZ4a/6uGYOmFCZcyrJKwSlmqrFt4bC5
fNkfk3wfDszIQvjmArYtk3m70TDZZgGPanOHyrnr/5lQIUth7NaECllKV1rF
1p3h38/EdOjZSYZ1tw12qPDk8cpzVPOtywTJN1K31VqjK9jGacG2PUOgB9N6
btKQWkE3c6hcEP8P+jjTM3OkIJmbFHRRHJ46msRnwohNPf7441/aa6+9pqfh
UJcYDXLY3OT4+2QM6tH4zTQa58dRqzOH+IOPhO6JEaqXP/e5z+G8GVaXkLdq
/rUUqwvO/j4ZYCyg4Ls9+QU3x1DAwDRTMVWIoW7T5x4O/vkXd9NXleWB61oF
j975a79LvfRAmMylwQXJAF8Ze12XlQ0mSWYpecLQLz700EPPieinJwN8T5wr
NnHond5qoa9UfawIqMT+36nNIo0zY6BJNVWcBeZofUI/54wZM+aq3XffXZ6Y
oc+5rPHGDj0r0EtilPuTZA8YIgyVamGa/Z60RW1lxZ4oGL6Kein63OwwT98e
JCVPlJLjAx16Fz4z9CwgZh10kMZy7sqm/yAl0e/1Cn3fbkMu/L1yzz33zNED
g3c2qBW9G17/bgM/cHRHAJEGFyfnZ05U+HULfW/oTw1SejUMNXnC0vyuDGMj
h569FXAhD/620Pf3FPpeF6eytGEloZFEfe61sJTYN6qeoafEoQ49WArC3FY3
F/y20PdLCn2/dhI+Kyshs4c+tzoM0LERvRUsuETbhg69G/eGHiEjbGnw3mRG
mTCk7/slHWFRJ08Y+geDVZrYBb1zET16N+7/rWB2lqaV9NfjCNGYiix1WCnJ
Ez31uYWjRo06P6K/UJ9zadMNHHo37n+TVAVBi1Dv12ZBAzp5onfSaEaDLN56
6XPXByndFPQhR++yjx69G/eGnlApQUdp8P40Xa7eaKONRPi9C2bH0F8enmNO
kyesiASv94o4W6/v0Du337ZQ4IgRBpUGH9RmQXNFDHv0LKheb33ugmOPPXZ5
HPcX63N07vCoeus59M7v+1USFutwMScPFTTf0BPrQSfkCUM/NUx6z8Zxf3FB
9Tx6txw29IQjCAVIg8u0WYkAx7VZZ4H7PgWzY5U1bm2wrkPvVjz/mqhlfSXN
PmwDKpBYS6mjwH1ffe4tQ51z72JNHr2zeobeGVNDT0fQIaJ67QXu19LnXgjk
HBzRX1IwO/0d+m+bxZMH/0UfnxK4JEAlDa5I4/666G40C9wb+mf23Xffw7qg
dyETQ6/rRFkeEufL0WN0JViH6j2aVI8OkfV2I6ko+Q6JOPfT5x4N1nJcRH9p
gft1HHpndA09M96d4S9pcKU2i3G9K1ZBtBWe61cwupcWuPfoXYD148LAe1ya
7SCoGiOMvIibYe3WVqxLAnXTumClmcejj+WxukXuxwVmFGt7zjTzMAGyah1F
uYBMbERplWjY7yOinNd2KF1lwy8LjD6hXtOLtqwjB0eGzLDNC17czC7YoGh+
HFIe2/jwMTPq6i8LDBo2dPeGgClSYdiuDZKdE+cOq5hzFQv9usf2UYG3J5MK
fz+IfQXpzP7doyJAOSHOkp+OyjG2Wpt9NfTU6mCG5F9r4Mot2j8dVVQ7aVC5
ajL43grovhuU5fRx48aVsFkFoctcr9U9tg/TuKknQGVMMgxPB8bunThx4uKg
FIbtii984Qt3xen2KwUxemxutjdspGpJmmbY2heHJpcec8wxV1500UXfDy74
+JNOOqmyMXdx+P7+OMUaSmflPUo3w/+zPo6io/DS4FOK8pXg33z9zDPPvCy8
ngko4GxNKN0w/3SUJG1IRmdcduDTfDfI+jJ7zZ86deqjo0ePLmGdqc242bTv
p2Il8Q6WDGsnaaMngrGOWGdecsklb4a1IjoUrWkJNVmfaXHlvAbUP5HP9tCj
Y94Jr/8Jr9+E1ziab08WBJuWyb9xf/hlcRiokFS1yWfjsfBv/Ml8mnMujpHm
XJxPJ+258Pqf7PVeeEXSov3LSXs9DPRrzjnnnJy0W8K/nwtS65Y0K6p1kdk+
jjS35lTSGpASyara+Gwmdaun3KeSY/JWoE0oaipd3wmu2lnHH3987heV6HKK
1OdTRbYwvHKRvRxe0RmBh9ppMdLwGJ8Pgzgn7b4JEybcEQf6Gkhz4dw1kPaO
fDYh7Y3wejC8QseNMKrcysyoothkxSGHHJJTdc155533elzlG1UI8MkuVLnU
+Rqoel8fx5Wiv8RBfFqbjVZNELcr+hfHjh3LaJMn+hRG36w0ZOpJ99PROxNu
6FcFqzcvmO8c/UPjx4+/PfbJGtC7vJdHP77KZhBDT/HQLTE0aOi/Fwb87Bkz
ZuTo55x11llPxZVeCf0VhVHU+1PRU780PdqNp5IT8GNisCDuUPSYvh9S7ML0
sgb0bv5cA/r39HEWB1Q55NxDzqth0hfEnYr+uTBVXx6ts63gXNzK0BcSbGtC
jy5Qr5Ojp6RndVgg5OiXnnbaaYsOP/zwHL3zHgw97Y8vo3dJO0NP6IQ5Jkef
zzkR/Y3nnnvuS2Ewfhp6N+949M65fa8w5T1dGPc9FP2PwphnvSpP9Ooevcst
rgH9z/Vx51o8o81SYLjsiCOOyNFHFyNH70InVyTjWieXPHq3mvh5oe/XgN71
vaF3vr6hd46VR++WUD8v9P0zhYm9R6HvDb2LV386ereANPTOlVwDesb925FM
Q++WclqLUfZhDksdXdcmedJcKtVIo3gupiQiac4g91TSvhus4dT4nEX03NLa
SCv6MGsgrRR0+KeCuXpWUT4VPKo5QUdy0py5MtJcqNUqWFp9mE+ni/DvnrEG
xuhyXqLR9eiwYcMWhDkrp2vFKaeccnPQ75wuF4sxkRV9mE8nzUWpjDScAYo9
c9KcovUsDDMj7c0wHPeMpOlniw9jVLlVoqfKlT4ZVW7WNaruDfZkSXBncqrm
TZs2bVWcdS1ozRT4ZnR5rux+9JUih4beeVzPFmyPoZ994YUXfi96XIbe+YGf
jt65zoberZFL6Hsp+h8feOCBcCtPWMLChW2v7N7wl6LGhp7K8ZPtFdGzkFgZ
0OXoXw3qQHYoR+98jk9HP77KJv13k6pQTSYNPleYdQ39k4HT2VFKnYVpz9Cz
Sp8UF/EevSt2fzf1VZ2CM/QxZJCjXxX6+dou6N28Y+jdGsCjd7n8NaBHF14f
OnRojv6hsASZFyeIzoLhN/SwN7KM3sVN3019Vc+6z2uzPw1OQFxiR80voXfG
3dC7uKmhX6tgKH/aPfo8nBHRL5k0adK90Rp1Fgyloc9LOrqgd5Uihh4HjSr3
HH3JTi+aPHny/TG2buhdyN4mfeajlV3WGoaeMrs6CWXoGQnzYsjN0N8RBtMD
of9z9HdOmTJl2ejRWtTe0b3ZcduqPHpXcGHo2VsxN4biny+M+56FpU5H9+Pe
ubsevfO2Db3zuAw9LtDLoQM+DX1p3LuKT4/e+Zsl9BPls/3s8Ho/vJhjPwmv
4GxdYoXeDVygHwXTkpPmVuAd3dsEVxnrSXNljEaayyIaacyxuQvwUXhV4xTl
S6Hfb5w6dWruKJZIc/bCVNZVEPd1pLnktpHmskxKWgNSvNdkdLmowRroenz/
/fc/tgtdbgo1uvoVRvI7JuF8BjeRrQivXGRvhVcUmbPwayDN1QMYaS5K7Elz
WdF3CoqmteVNSHs7vJ4Jr2DaJo0r6JhR9fUZM2a8EjmzIgmiPN86YdiwSZO2
3rpqk8/GiBHDhtX26PKC3nlKXVlzidIXFZVzpzsLI1KLZppnnz127LJlO+30
6w/79v3k/Z49G2+80bcvHnWPHtOnm92xLAwRmceiyq3lqHID9Z2ClTSqbg/d
8lAYaN1SpZVE7fMC4Pvvr7MOlN075YADKOuxVGDz4OAaP/fc2mvz3muvXXYp
pd48dW7htgbqGCAvhl8z6nrlU8ikSccd98orgwe/+OKWW/I338UsoAVzZ4Z2
n3yyv25zuqRg4IwyK4N0c6sFZJ2nb5QRSnojzLlEs4yy9kjZJYHsjz5aa61v
fGOffWIplMvMGmkrVqy7LuTJvyxX7VJdnjRXIGSkuWn/JUXl3H2NcjVZhDz0
pZEj3357ww1RMKPKefyWVacH33mnV6+jj95++1g94WxIP0eVmw+NKhcNfqkw
IVnsbdV5xx770Q/WXReKYtGWq5qyOoPZswcNevBBO2bK6klcvtOT5ObInyQO
au/wW4rK+UcWCr7qqvPO++SD/v0XTzzwwGgkXEzGSELnP/qoowNpxQIflw5a
25HkCnzWQJIbjg0lDJJ+/vPNNnv5gkMOYbo2qlyNjx7+2Hj55bXWmjbNFiXn
F4ahp8rtVTGqXCTAqMqTVJGqt9/eZpt77z3uuBjEMqqcK9sjDUF0PafK5R6M
KiuCd8tao4o9w+NjpE6p6pkbVF7Lwirvu9/deWf+jimQqFUo+mu26I7FM3Th
b37TroPxvIJFXceR5XY9GVkuyKxk9cgnxJtvnjLl448HDJg166KL+PdPg3RM
UG65qSuJxp//3NaGQsm/rF7Snd7hKXKxR6OIbagEegXJy/LZiZ7fNH3q1F/8
YuON58xRgcXJMOo4wZNro44bOViBceNsgT896XhdxtXfkeNCXkaO84yVHInE
/mzZZps9/nhAb6m0mPbUKVCfc+UcRhKW/KKLLHk8vSAhI8l22Ljynh8XjJOR
dPlFo0b959udnWLNjaQHJk6cWKeqjKRYSJxIaps3L6hAQTbrOkKcO2OEoOPU
luSEvPXWwIGPzNhhh0cOTrJxGVAz3G7R0rswzM4tyMaT5Kzkj5P6s/k71+cb
b9x//zfe2HTTfMaDsJeDxqDEJhq3kChRJDm5NtmKIxMG8z3RKnWFNevKd4yx
yp6LM7E08EVt1G168iy5sIBpIBuaJ8dYi07ePSdMGDfu3XfXW4+Jmm/QSCQO
W29ZKUXscheSVr46c74QGbMwA09ObWjDpge48D3uDG2rI6NrQr5jpWzPNdmP
j8+NSFDVk5VsZ7XXc1y6javGpVtNKJc97rln991XrdLCcV504Kuc/ximPDJe
NvJd+NYYvOayTTd9/cG+fWkWfcDFiDvCLdPdgZ2DSSYL7CDVDcRIT0iDoT7t
bP28t3zg0DTPxS1fUMrw6PCjpob5Jyo+WF8LEPinsYtcyaEaiva4n/0nZ6+9
9ruPhq6mA9ZXXxaLTR0GORq0Tjaz2+40d7CbJ9s5+j9O3VTnWI3sXO6RxO9x
G0EwHvkYdmk+nSo7kBo1RU/ZLtoFgU6c2420MbYdUHIIZ3B4dLIk9UmgGzii
3aq4NMqN6E8+6dkTOedErw5YXgu9K71uRLuVr/pCnVhU6nIphIpDjfUChFcD
tUXagCEYg0E7i9Otq4zy9bq33W6+tRhftE8Z5YKNbMtMK0xvFMy3ubvs4mA6
kCJze+7OO4cPf+YZK6nWNJ5sSIVBxiq9ZAeCuO0MGzoGXIS4pObGAMuqaImM
gR6o6/1WZx6rq/klescdqSd5DL3QmbyNBdogbZiNCw/HrVq2MdYd9DfA0eti
uiX9tog2tEKzDiall+E0y4q6cr923uGHHz4rWDtb8RBseCTuprTAAIJG4PIv
o5w+fiIIWmaBwxXW7dzwlOPH4niLRfxxsj+ogCivUY6aoC455WSX4qqKVywo
j46slZJAycrQeVLm3llobjOzRWFiZqOTbNMdXdAR23Wxfvem0Nnyp7VZhiZD
NKNcUL0aD9+1qtxlQVd+EBNV9h3cPR76RsqGLNfC7DdzppUuDdbv4tZeIT07
8aU+D3CgI92dTWakuwo3I/2++3bZZcUK87fakxWvK3qNzOlBSj8Mc7tY47aE
/tFBgwZJJaflzn796z59WDDKvzZPpC+z66Sqg5MC1NusNnakl4yiW409XRiZ
cdkVtLleBijpcl/JGVVKhbtthUY2Psd//mdHB/+vttTvYOvReLqIbuj3J7Ju
nBtGn4V8O3HZtSbMDdA4ysKoq0N9SnbHKrMD8YTPuNmpV6GZIWlg1vtf7MQc
d87xJo5iFxwwQbvqOiuZRT/Qk1zQBAyWx00pRrGrgAqvJ4JElgd9sdIRAi6E
g+RfWyUjWKvaAUmt6i1wmzqKXUVDSTWMYjdfGsWL7eChXDUYISeZ1ULtHgpq
9yAbINFoEzTOAk6D/EsnncaiYcOGTY+Tx/6J7PqQcHOBN+x+snRVMFp03bj2
2pEjWWPkZDN4cCpk4Km56xFVg1e0UQgPH4a1hyk1TdFkTrnbKrFfmqfqI9E3
c5SXVMSZEaMcZwKn4lPMSHseqopB0NqchM8nwmr01GyjGGMyhg6rbfQ75vkL
41HB2Sn1t0dvZbBjoaQzrgboiYLlMhbcvhZj4Y79Ntjg4hFVNWyYXqaUb7vA
U38szPZSrhYzWPnksG2hU42F+fHA/IyFAd3rj5uHdFdC4/zzjzgClzxnwXW4
SZooBuFguRuyTT4bH73R3j5mjDNBjfHBObs/zLTCjrXnHDr1PPwGmL0V1m0e
3tyxU1IqF8U1dpyrbuixTTdFN9/YIfgIE8qIfhKW381ekZ0zjR1RsGZhdG9X
6J3SdQtbOHZKCuaKS42dgn/TyP2byM7vftdo1Ow0lZ3YM11755GgbKdmCwCn
xEP1u/lhFXNpnBOHp/EyP96GYYnsjbpXNhe3WlmYMmPChmNQYswrDusw2Gp2
2pUdUjk1O2rUOsba8I+1VbABOzvYfGc+tzy8R5KyXFYAF1s5Lko65kI6xsVV
Vx10EGH3vFOcgTFh//rD9vauneJSUepDt2PEHt97772ZQeJCpdJdog1WdfWm
+d2SJcL2ih+0tWPB6ZVFElz+z/Z/uem5WRj15iQ/eVa/fl2HyZIluqKQJ+xo
JvKATIJvvzZwIPpq1BNlqBdMVkiGz1JfLGPnddihqC57aR3g6sUeLUxSzW47
oDPGfk8dX1UffE+5IMR/660bbxwjWlD+5983m3cee9BBUeio7e1xFrJDk+Jt
TzKAt3Vku7PES3Fk2/PnPIxGt0Jvz0N8MUVke4PbGQXYqSvCQGZBI4Et9cRl
+YABkbZ0LMrpNKyz4t1W2znK3QnGRrlbID5SmI9KAjfK8xKXeNyDUd6BxSXt
w79Yhq0KjgelEkZ83KomD++YVIW99bJItKsw7Zze0jYLJ3Pb5unco0ZBTkZ5
7ofGJErcj83Uh9pglaKFwnWiyOO20aNHR6OJR1kv/XbwaiPO3Q6OhdJYdcI3
Ftxyq1GYPY2FPH8W07GRBfLMzN7Z8w1UiPICyaqYk0HcpBarHkMkx67lKmQn
PG1SGLOljliuLePF483Lv9oKk6aZdFd0UekBW0jJGHHpgnj0WfhkY8fto0aN
wu0wt9Wt67RTJHrOyWKiUTs7Rkqj2GUWjBGcJZymnJFCd3TErSXBVmw7oqpG
hOXkMJaVdj0bEwG5Kkk3WlrtstAjL4RhIbvUtKfLdTfqSUne5dro2Gxjwgjk
Mgzj5YrDHJNuwL9T6C3dG975l780ZLI+2V6YKU7DkJ6y0M8/9uu34dNVtZzp
IXiEfDY+7Oz8bnA7NDRlGR6XM9BMXQfDhx2GnKYSNzZRxtOQQhkpm4moLVzn
jbLOKPUVj4xc7I6dC7dp9zbCzecPJa8EhvMedau1ZoJ9p61tnrDaVIaDKX5j
N3vFNCT2OXorNs00SXU8FnT0zqCjHwfumjC6acFN2TJxFi8m5dqL3R1nLtxh
3YiHWx/cowccdER/KMoyrsAt2Ii6/6Gq/l1YaVeGYu/JE5YaJ2VI6lD+ZQOR
NOcFYYT98+LFiyvrNqqIunYbRTP14SY2bIs3t+6RLFc3KWN3JItxmPtivOLW
7mxH/6rOzs2Etw7lMPTMpK4cYhvrWhyzkqsXLlz472+++eZPV65ceeGXvvQl
6y2Sn1fFsa5DT05Bjfev1tfZftYxVCrfcQeCPKDjPV8rRbMSI9vxpOj3quq9
yE0wIyNqbuw5V8qg4645zYwkOVjq/a17XIRDx5bMhVgqolP1ZbR2Rqddl+Ac
fuPG7dPQAz7aoztLl5EoEIfAzHaeXB0f/jMFDN7lkq7dg8t5/fU29uMRzeGT
M0FIASyMxQFa5ODX4ANTH8UrpLEXsvdqH8eVmwLeKRgM4UoC8jpVTok4lKWe
LjlrbAW/83tw9n7otJ83m5dGK+EmAF39iFk83zqK+iJy+Hh4mubyK78uZ2OT
JY1XdeOj7uc4c4d5W72iO27iPsUAa88s2WorPaden3N1y+2F/rfn6KC4NqhU
hh0USFGEwAlh+Bqcc4D29S84fTo0xNsla8ykv9oOUWaloxfUxXtB3P0J1lXu
dCFj6OGHdtzxn27s10/OWLHh4TKxdq6iq0HO6qZqv09Heie+xiWmeZQ544NT
1oonW29/DfifjQe+rWNwgSuq7ykZet6OuR4fp+fPO9bcHG195c5tsqODfjB3
gw1W3WrzqpHssm/GmisqsecWLhwwgFWQ/EsHf0+ccqpzcWvfsv7CWsAemU1h
0RIlLqSiU1yvrqdgQ0s80ZtNB8ONJL3VMV4k4yZp49YdkmXnNHFu03dCR1Kd
ES0IqnNOePbR0HmM5CYSrXqgMQQMeS8I7w+CMCkaYDMH6RUrT1gxbMCA3/1z
s/nVbShZkJscmvFkb96ozOap1yeaeEjwRhFhe6o2dKDqEUXF+80gmvnBL4nN
IZamYOgVJbOmM8/pwcOdlNyEb/sT3ElmepZX33iS2YIgpQe/sd56ECCDiWOT
w/eIgEH3k6APjwcKJJaqlPWIt1DE97OB/1nmF7JkspKHt4LHvuLK7bfXayTk
sze2kNUPb/ISoMGtXGKHwk+zTE52VmFvbN04e0fZzg6/v8WCvUqHzMtI1Z5t
RuHEs/D1ctl495EryflpwSbocW/94jUR0L40uDC/CNR+Jwji4yt69Xr3m+us
Azf8RtAkTmgxTRudWlKg8+KuC1uUufSKrpX74FhTr057REO6uzID2wOtOqL8
BtFTHJclVXDmYYFg7h05fG7o0KFskabqll1TbwSTvyg4ZUzJ0dmoY9d5gZbN
U5hw5qKcU5ds1NMlenW9fIOjxfJrJ+AQk6Y2x29/Ng63LPTjuwWTcJsuUugb
vDF2XLIxJ3J4y+TJk63UqmPR6aefHmuLecU+KtWmqh3roIaPMqgIEgsfNN7X
Ge+3gD9KfRhW+IOSqbSjE91lQV90zLnuy/bY1SP5Fo+GPTdvh9+ptF0+evRo
tpbGjfaU/H3r0EMPlVi7LVzcYqZR6FbjktDpalsusExgWQvsxl626CujWVIH
dih/PKPQnT9nXA7pvgsRyorAiQxFPQi2k79BQ0EjG+nZz80Oo5kXX3xxLF7n
uJiXwlDLqj878lBaPjxVrE2Wliyio55q/rAD9SSyjR/FhHCtXU1ix/65sxT/
wfFT6jWXytajfDuwkwwAKlmpr6XHKEi2nWpsUX4q0DUzqJeNuQ5WL6xiIiux
s+KV7OHzZws320ycTk1YyoUn6EbXC08stOAqcY2PrQr9YnzkZ6vaichyVtWy
YI6+H4wKmw2pyq36KhNPhz77ajzcIpWoqLLVdtJkOn78ttv+8od9+tzMBSqa
jOv2ApUeiZ763JfxjnzXDT/Tx+NBtNlx1nIHDUXznE5B5e7CU045BZVaSzl4
PGjSdXG/oqbc2vOlcuyCShFJjotF/bjw/6fDlC/ZE43RuXtUMOZySnvJGzNO
ti50hHEyPlhqOJGG9VzwdpQpDg62AEg969pWz/rFL35xQQzcWPITbZo7c8QI
6QMlv4OVBsv3mBAiCbQ8dKqGtToQHGgwaIwKTC3dYZFOd37wBMcE3SH3sOOp
2pZ4t6a9XD7lIhCO20GTGA39lX5Onbo7Lhgt/8yC8f0VG2+MVxBDDdAv4h9H
BanNB4GHO+JJPTukYUFUmtguCedpofc5IXKiEm03ZbqVgxF9eSCuPu70YpU8
YmGQoUMQzllB1bpK+V2TJk16oMtBC7OPP/LIf4u7LPoUKLeQxM3BAH09cq1i
kou37gsj4r3gEE6KAZEvObLd5aN2KsbVgbT69OQLUlOQTWEvFbCiMOsp2bd9
+ctfXhlsoMwXvYzN8PczwQKfFQYqpMpy1JwcNh6wS0nIN51mqlwSlEtWAzua
FQ84qSMYbUErKrcQc/9M7+0+VLfEztio5/YZiY35ock3gzrVbKyfpL/KYr12
xkg79uU6K1ZmicBC1MKqTVhiwcNeCtZyZmfI/BCSE04sus4cszI4O1hdC/W3
o90oDGsyuBqQTXd2m6yLVr2XPJZ6IJyTuLozoIQrNoVcecEFF1QbKFccpPRA
l4NbWPs/Fzc52HiEE0YzltUik010Dq7pLJZ1WW6LvkLdYqL6sjD9PLL33nsf
tlsdtqg5ZD1ycey3kx2HbrvL+2mo15b3LB3qRGkpSb5/7NixcLehMvdI6KqF
8TQHO5NnUZjrz41nSduwRtN+8IPevWWHXBufbfAU1hmBu7rUoE9SRUYTdiAW
etwUWIVlWDcPtsn8Eq8OYVaxLehDC5b4/WSJ8adlSvyycsasxbTOgTto44A0
l9w61rLnxpnzKmXcNHQjIv20m73iWccXmjWWScXKyNDse4JDiEYbt/XNI3go
TCrjHBfufuYPtGV8NSZ2cVCmKBf0FRrIHnom9QHJMXk2MMisGA8+c1tPrSKH
qYMZUKZzK9UjNoAFoVTCpn25n4TZjzUM2nS8Ert9wQR/kGZ9vCnxDvVmH6kG
wo6hTNVGyRt8IUx64kIZpe68Ql1BNVleIj62wkOhVRMz2HF7WcWTd9ddhLLe
eiDeLnOcI9bduvkLfRzXj1UlTnvct4p+sBLB30A/jF45FCcM6wsyeuneF6xy
2W4JaEIMREEcRFp9OKpEkoSIA0ZJg5+SHkAZMDxogO19tAu+uafillhH+osE
wToDQ4sSjE/6gJ19JYifpcXANDU/H4Oqdnwg6bXL4n4QTfu2QxKkIWWCsGJL
11VjB+nE++CI0UAKxI62l3uF4v016Ict5O0idBeM++dkThmSoh9fNOEEPfhW
cP8wpbh3KvImKoKTim+BqPulUflwNKY63DuR+vnGA/uZ7wsdirYI373T7IaG
sW+EgA+8EAIcm40L0rzoj5wIfZRjxoXzPkwjs9afk5W8G4MN+FbQMRykyrqB
8wi/HZYQuNvx4FG3u8NStFi7yaZAxI7bCE5pXEq2uxCTonssVtXBDlEeYcMR
A4ftMCjj5jZvBF3B2HPdTbxu0O61JxdFxk0e/LCsWnb9G8Zzefjq0aC+9M7A
NM89bechRZ6+EQbI7KhaUnDRJg5gA8JmWh9V9l08ZkcetttauA6HKJrMXJYE
1/Su32SbHfJaK9Nxph5BmV4IankHXsXApEkckHSnaZKRC+mroyZNTF1aa8AY
R4G7z/2XXgM4lNuu9ytqAPGL14K0CJrbwPI7Mk9NfUCXYVqxWIcpBbr6aMR6
YHn2lwlitl06x8J4bLnXbExhxtgfRiwFr8RIuSYIbU7c8XK6sjaDHYKxUma0
owPK2NYi9Y//kqYKogv4SsxUYwrdkQ3sZWEpxyFj2NP1kpf0lB1GUE1J8r03
zIdyK8goJWIXJcLtKbMTwnFlyE+RD6m0PrJxXVgF1t1h1pzjRuiOuUEUlW63
aSwNLs+U6MeckUbQXOuOysJcuxb64l9TX5COIdFE5EWXNY0FgRpmYYmG6Ezd
JObD8ZIcegLjGxQ6YlrqCO6zkZW0XT+5W2FQ/CrZZIJZhD+RvoajGv8YWn0g
9LqgV3euyXEQHPrDflNM/oZpXDweOyK7VOfa0OWypDnIEYGC0EcSFdZbAiQA
T8SQDB5aeYDNmMGsoJV3M4tumCQBEZxLwV7FjQpui123c0rw724Lc7gQYdkW
zag3qAS5ISZr9KIGuZ8QlST4BRGaU5e6oSeCQohKbpDGBRumkQREWKERZUFX
RINmOzdlzgsrIxkXdr2mJr0bTINTYxZWL6uQC3YgAp3AEbKdZgyRZ7jPkoiM
LoiaHO8AEYwLZirL7boNybb5GM/lrlh9N9IRgXCgT/Yu6G0p7dCMaYcIVvHi
ENuNMBhPRijHq9gKs8n4QDvZzAclmySH5xk7kyjuF8efoopRTIMl94anESo1
mhjM/0iUEF9hcU+8hbhLtXuyFTgwskTX4EI7Uyf7bDkRkrEinoltRloSvpsW
R+wFacRSDY7vFO8V/WxhxOodPlKmJEc1BOtJHIt4Fh7gMMUgJ+WGFmWBp+Ga
DgYNkx+HAXLIMDuAcfptn9FVQWzzu9xxARYcOSm53V+J2lOJcvuM9O6mDkRK
LJCY4GrblIG+7JCIeirogxClobAOTCqnXNBpDCT2RzOMNy+MaDsFxJU17+uI
YvplL5r4ynr/l9h6UoWElgi8EoCFElk1bpM6sBZXH6UMpeZoLShjcyszD52Y
LZEesaEez3Rx10BZBHSvgkr9WSmLOVn6kHwMw53MheS5tkxK/mKYmST+oEcb
daJa7CXHM+DcOZS9zgFvXVCxy1Nv1hsOjLy9lTz2B10Vt8boPX096F1SQxgB
STTYzeeQd1a9ohGr9HJYzYnya2KgJ6YR+UElhoFBgMrRu3gSEjmxkj+CAzfH
mWNmGgz17g67Klb/J8vlCbEwSm907ANj8d55qI33teNaRKHm97VrOcNaaMDz
oQ9uOeOMM6AWPaTHkSsH6nG2AnYV+UI5BoUBs30ybw9ELyw7ZeyuuKFmb0c4
s+BiK1bMb06VL6KsUVQ4IJr7puRGde93TAlKurNNU5+EOuBIrn9tynd9YoqP
92PBoP9TsFbEgziXFDZQ4nnWCaxHhiUP5tFoH7+aLHW9pcmqb1RpvEOZXSQs
GhbTt3HHfUwBS5K3TT6lNh82Me2Vpn87Yo5W4ulhBl6w7777oifDC07nVWlG
q3cpqwkwg9BgGM5199TqsQQNLjuqq06s7Xn5yvLqNGVfG9u2UiqtaZHwwQx3
a7ttKXNhXGt7VlCmOjZwTfJJ6ov+7BIfta6NG4NZPcVddK11E/4Sl+GFFeV1
+hyE1SVAe7i2nfeb3W/vb2YanizvihgjuL5g36wicESyvAuD5ZXJO7vJ3t/m
ZG0ziFba8UXVDWnI3xtHjrWts7KE+etbN6VtuyPN3ae1R7LLsv+HaNBN+hyG
dm7sS7tCWH0PudGTKzWTKltExuURjG5CYmdGPbm5IG+r/1TnqsEkdX3c/ixt
r52GFrYpl/fcoINXRjN4qz7n8je6ejDvscHljRfGskZpe60k7wdC27jxsW03
tX8jWap627a1re5xg5zFpFhZn1227seOyZtJh4Cn6OBcfc4d7jTMtc13DCt0
PL8j3V9KuEfSb3Z1yOz/TX0OR2hilLdVUB/cXdu9krzZ/C3eY6nt2wvytiuN
D+lu7Fh5gqs63L3gtcwvyNva1oWZOL6sDlLbVtFGATCFwF3bXhZN9R36nLsF
z8rjR3XXdl6ZGotZS20v1Odc9Z1tINBFtayf2EmTxnyzMC53S2P+EQsdVXfp
cy7Tam0f1l3bbWks1/txSm0v1ucIEh0ZiwNs78bhqe0lzlbZ9dSuEn3XQtt3
J1v1QOxLa3tMgW67cNpd6VNqdok+53b1280AR/hmpRv/mnqwvi1s19SDdbP/
qM+5ewZs/82RqdlaO+x+ZHfX6y5eMaTZpanz6snXmh1bULo/J5tQb3m2Zoly
1vp2vz7nijpsn9NRqdkFcXz/VzJ39cb4XdLQJsUkzT6QnqtPUrJmv5CGdm02
/ljosmEFi/GQPucKbm0z2dG+WSb7eI+uOzpIb6CWhFjd7LICtdbsMalZfAeM
f7xP1u3asWa/dsQRRyyKzT5coNY27mmyQeYT5qpsyeuHw07JVZobphIx948U
qLVmNeEiUyBLDpm6/71A7Y5p1iaOhRdWrSxQazslT9BmMdJ4BDLKflugdoc0
ypYFQ0R4rHosUVtv8LFmNZcljgZOjOjtbwrU7pCGA/6LxINXFfTWDKYm+sQ3
Qm9ZtsTrIB212YLi8ZgafSIN8rqSy5o9KWkCqyE8kBhHdNRun2b/B6ImrC5Q
a3doaxGceLdz4x3zdp2cO/VWG2/SZYvC0mO+acNTBYKtZU1Wi0++qE+fPrIC
+ihZkbq8fbu0GmFpPR1leLZAsO1l1tIiWY2wsMJliQkOt/l32+SpsOlPVuzP
6XPuhH5rdrx8ylYzPKFrbHfkB8n217ukrWWSvmRwxFd+oWB1bS+5VlTIPerE
AiRKYRlft5d866S9DArJwr6UnqunemtW60skxECFhiy7rAjA7bjfKlldtFcG
xbcSV/UphdbsqSoHrC7OHOk+5JBVn9XHEVjLLOWpaBAF/rY+54rBBuv/JinB
J5scRB+sSM8d1KAxkCZyIAdJlABNeyXpel1Gby2fnkRxZxAFM2qsxHUn+av7
IPe6kk6W1OBryT7UZ25as5OTKJigCYkhsdIW4c3T6CBhSZ0+avG6PurOGLVj
K85Qmgn84VRIhaOdxeJONhmcuo/zGsRUvpkoqE8xsWanarMYn/nWfdWPEov1
dYaDfd/JFP89fc4dHWvHg0xLfYeKC7U/KBgfDXU2oTbvu+/ro+6cZmv5TBXy
0IxgCP1+YURb6Im0HTRL9/2gMKKt5bMTzYQtZIB8N4miPtRpE6WZiZ4oPhlB
aDapOU/Yzng5R2keZGOEQD6K8UbB5xnoFYMsLpbIetodIGyNT1eyR5liSGXf
64VhslGaPp7kNiYUI9vPW48Ra3aG0oxikOgma4cleq1gMTZMA5CAIsF8xna2
36W+vssaP18bJ6jBLEA1FiL+dkHalrShlOP5MOsRrYwloU7UdgDQBdoyOs6e
PRwOuHtZIdyd1Rrfk1oLKsexn0x+PyvonjV+UZIJ4mIhyMxkRtedsaaZxyY+
AHWa95tM3ktdXi99rfFLtHGcFlY7rEzEdy94AuskgaMkJHnJ735QcF2s8cu0
cZY6GH8S0Tz2bOqi+mA+a/z64Bli9C4z9bYCB7f4tWOiLtfGsXvk4inwIUz3
TGG8a5FDkxgpFQKcbozd+6ggc2t8pjaOgaHoa4U1/lRB5n1948TTsFMfF2Ru
jV/R2jjkPlmQee+kLXnj/1qQuTX+1dQ42oK20vjjyQOqzzvulRonghnLaH5T
kLkdhjZbGx+SNY5yryzMNT1aG8cl/W3BpFjjX0uN32CNI/MVBSekI8k8b1wT
c/5CCWv8WvlsH2Iy58Sn8+xYyWq5grnTe22HORaRXBQ1PAhfvOrfpWFdX5th
aL6uaBir+Kps4p0dzyTUZVnP/MgRXniU80ePHv1ysDNopqwM/1CYMAzFDYoC
vcfp5Gf2i8i6QBeUndGdYLLgCEt2GF4aTJh42bY4dVee2DF9N2sPEPrBGKAB
sstfIKQS8t211177CatZIpNe/aVgFK21W5VQdIhSN06NwMUQj0qDFe38zXfP
BSX+UZj9xP21EIJzzqzJOfLZAe94KpyvwEQfi+Q15NQZ84/8xCM8SpiZ2aT6
b23dXbdjrX9TCR5qOvJ4TvA8TzC/XWybLytVNH/EsDHh7jQyNLe3MoH40Fr6
+1aBlHxl3LsamQBEIzD+PJ6/pvFXe1qG7A5FRt8wJJZ1QXa9R0aPr8z2lari
+XzqXwsTrCG7a83Iru4eGUrSS5G57bPWW/kBbfE0zLvXjGzWmpGVboczZG51
ZMjubUXGESJ0N2ZAb1qQqhEUCzpyZHYPhDuA/m8FK2vI7pNPORYFncLGMFoZ
ELAixmS6qjrsoWPsnKRpuh+1lrFoJyi7k8r/VhhchvZB+ZTJhJ8gviufdPpU
aaFHTKfl/KLpI01P7fBJd/eFIXfb2wz5w4o853m19SjjDR4Ze1pQ1gMRYJmg
i+5CPHMy1ksnPq0B+UpFDiTDDDYYcmgdgiIPDqsnKXJ6GS1D5sgVIWHgoBu9
2bR75O5iUUP+hHyKeiN2dgPDCpoJONyDDQLSSby9ovj5iUfyDcVo6EhDYHtk
3Wm3ayDmGSVmmEmCliAG+cIqFhQwPAq6Q0/xlJhUvZvW9u7GLqEJmqJJ2xrj
Nl3/rTBnGzEvymcH/8Jnh1a8DrqG1ttAJUlbyeZqdqpvFEzMuaMiPIgmISAm
NnprfJW2ysZyZmfB/1aw4EbYK/p4R5V23FZtNCjCavu2I65q5puHEZUYbfEM
q77xbHjeSDf2YxTfDbb1+RjDIwPbivX4vr4Vz8h1m2ON3O8LfcMikVUz3wSO
PMQMiinSibFOia+jKfGoYPKEcFY14/7s+JYGq3aaHm+dNTzKZq/u6XX+sdH7
rjwtkI0TrS0hfB35bER1GlbzpUaufln+1B1JrQU65QtOL+ueFk51i3n7UnvF
W0mtPXeZZaE9t73fCoharxK1xpwr+HcSV7z/09pz+a1Ce87BUJiWSzv/Trn9
vf3+aTCloZjDoBfXxUKhvxOmZIvcTZAFGHdczRqMaQ5Tmg2cD/p3wpSmrxym
NNn//8Bgi4avAcZd8rYGX0ZgrDal5Gy5hUMBxh2hZjDubvQCjNu1lCUk6ysW
N2qFcffdGIxb4BVgnL4ZjIuBF2BcnxqMWwwXYNwFnH8t9E8OY/U08YRBedqy
p+4ukQKM26lnCzum0+e7rEBzGGfD/lIwLQUYdP7KWB9iMC4YWIBxDsxfCuO0
AON2RP+lMBYKMO66kdLiNocZoTBuefQnhXHLo41bYdz1Gn8qyKAA43TnTwXd
KcA4fv5U4CeHOaDAzx8L/GzSCuPOxzcY54sUYFz//LHQPzmMzelOry1z7VIB
BRhnD0qRlU1bYZwr+IfCXFKAcfz8ocBPDmN1WG61+/sCPwUYNzVbMMrNp4Na
YdzGDoNxCZ8cxuq4nO39JNmDOgO1WSuM0+tPCnpQgHHzj8G4BWoOY7VaTkd/
m2x8nd4d3Arj5vrfFux1AcZdtmEwLkCQw1itl7OJFk11vuDmrTBuD7HB5Kfp
l2Dc7svfFHQnh7FaMbe99FcF3dmiFcbt0vtVoX9ymEMLeP6lMC9s2Qrj8PxL
AU8OM6og618W8AxphXHj9JcFPDmM1ZM5G/JhAc9WrTAOz4cFPDnMYQU8vyjo
9datMA6PwbiahBzm8AKe9wt4tmmFcXjeL+DJYcYU8Pw8ya2uztq2FcbpgcE4
3zKHsSI0d3aCZeoourgizqfbtcK4usl3k32r64WGZjBWlebm4J8WZJDDWMmZ
q1nNzlqtdWf7Vhjwk1iVBMVPCnYnhzmqIAPL3bq6vR0yGCswcz7fj5Ie1OWg
O7bCuLnEYFxJYw5jFWduF5elw128eqdWGLcr4S2FcaXYO2cwVoLmyv8tPe9K
w3OYYxXGbXOz3LsrnR6WwRxX4OeNxE9dyr1LK4zj5/UCP7tmMMcX+PlOgZ8c
5oQCP68W+NktgxlX4OfbBX52b4Vx/Lxc4GePDObExA91VqLXL6UxNyfe1jE8
g7FiMrc/6UWFcRutchg9mkNq57mQQZ5+XmHcaT2fzWCsuszVxD9bkNueGYzV
jdHezDjmni7Ynb0yGCsKc3s0rOjNbenYO4MZrzDOhjxRsKOfa4VxxaNWCuh2
zuyTwejpRb489rFkD+qNQ/tmMF9KMIujHX004al3Se2XwUz0eIQ2q9V0tx/s
n8GcmmDqqtXlSQZ3xo2jIzKY05IMLo2FiA8lWS+I224PyGBOT/1T96mV/oKT
/S8i65EZzOSkB7XuWBUyZbD1VpoDM5gpCjPadFT0+l6FYaPblGivP5/BWCUW
en133G9ydxo/4JaLtA7OYKzMCnky5qR/FqU+pfhVxtAhGcyXFYZxiozOCj5c
LLfHT2ReEHN4aAajdVdS/nJeVoK2QMHYFTEnnj8wKgM7S1FRXFaTNy9TuUCe
7KY4LIM5u0DebQpDuIYyEyHv8AzmnDJ5ttmEWi6iPKINYzKw6QXybk7kMQPK
kShHZjAzPHmidTd68mQGH5vBnJfIQ/Ni5ZrumZJzXFkYRBKPyiCt7ChWpTEH
w5jmWeVQHEpuuNEOC3l0Bqk1RXJixcOxFlE3rjWhkxlnsp2WcEwGdrEijLWG
kIpYNNcqlyHg30IqEjgug7w0kYrLgZsix7irbCAT2YhBOiED0yIf8WzwhvD0
7Vp72QNORGuyHeY/LgP7imJD41ggxCIY3c4sm/mhE6WEzpMyyFkKmVe4ADk9
4UQtYQFaT8kgr/SQsWTorCRVIMGJgv5DBnmV/NnetdhIzIKemSC77qEF0QI+
IQO/2oM/YqUtAj5JwWkOcAiDwIkZ+DUKHktWAOcsdckyTxDw9niOdjwDSwz0
aVkb18mfMsJPt/Q4zIIWWL3lt53eI4YXT64Wm3161oyVtgwzTij65XmYr47X
NmRkBBtJDJiuFVN8RtbGTUoKmjvSuoLY6dV2I+BR0kxH1wPPY3XFl7OWbk1M
MQaWWi48Fhgcpi3ldRr5GeNnZS3NlT8loR+rB0AMH8hE1u96vaDEtuAXesGW
V2NU52ZtzpM/e4AK4pHFSqOQdmmDtj4nzfaIu5Bjtp4uQDQwBlnnZS1ryYik
4o8wCcZUfKwDoDndvt0ZdwfT3wSK0dC6AuGirNlF8mevIZWvakAQwEJZzHMj
EF0xyvqHr2KeFFMYU7lzsnPCL8swLZE/+wywn06sUiJ7pWUrwYhq0WQ8iFrn
+z5oWEweEy2KjK02jDHDOdyYnJVhfjApMMo3vkpHNVdN8DY1yYr1qZOsnZpk
zQ+/tiRrr+ezBGvM7seM8KHGHpp+dUbDY3w0NSu6T5VlTzvls5GfEx2zp3Uh
GK8bssYy1/1/0RfVZ/4fPBMbdQ==\
\>"]]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.609778401970359*^9, 3.609778402004361*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.6097776155213766`*^9, 3.6097776156093817`*^9}}]
},
WindowSize->{1904, 945},
WindowMargins->{{103, Automatic}, {Automatic, 149}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 934, 24, 31, "Input"],
Cell[1494, 46, 1714, 46, 52, "Input"],
Cell[3211, 94, 702, 15, 72, "Input"],
Cell[CellGroupData[{
Cell[3938, 113, 2015, 45, 92, "Input"],
Cell[5956, 160, 26868, 446, 611, 1301, 25, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]],
Cell[32839, 609, 92, 1, 31, InheritFromParent],
Cell[32934, 612, 96, 1, 31, InheritFromParent]
}
]
*)

(* End of internal cache information *)
